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We use supersymmetric (SUSY) methods to study the derealization transition at zero energy 
in a one- dimensional tight-binding model of spinless fermions with particle-hole symmetric disor- 
der. Like the McCoy- Wu random transverse-field Ising model to which it is related, the fermionic 
problem displays two different correlation lengths for typical and mean correlations. Using the 
SUSY technique, mean correlators are obtained as quantum mechanical expectation values for an 
U(2| 1,1) "superspin". In the scaling limit, this quantum mechanics is closely related to a + 1- 
dimensional Liouville theory, allowing an interpretation of the results in terms of simple properties 
of the zero-energy wavef unctions. Our primary results are the exact two-parameter scaling functions 
for the mean single-particle Green's functions. We also show how the Liouville quantum mechanics 
approach can be extended to obtain the full set of multifractal scaling exponents r(q), y(q) at crit- 
icality. A thorough understanding of the unusual features of the present theory may be useful in 
applying SUSY to other derealization transitions. 

PACS: 



I. INTRODUCTION 

Derealization transitions control the physical behavior 
of a number of electronic systems, including dirty semi- 
conductors, metals and twa,|dimcnsional electron gases 
in the quantum Hall regime.Ela In three dimensions, such 
critical points occur at the boundary between a diffusive 
metal and a localized insulating phase. In two or fewer 
dimensions, however, the metallic state is generally un- 
stable ("weak localization"), so that derealization tran- 
sitions are typically isolated conducting points separat- 
ing two localized phases. The prototypical example is the 
transition between plateaus in the integer quantum Hall 
effect (IQHE): Within a model of non- interacting elec- 
trons, the localization length in each disorder-broadened 
Landau band only diverges at one isolated energy.13 In- 
tense experimental effort has focused on the quantum 
Hall plateau transitions, and has led to an unprecedented 
characterization of the universal scaling behavior. In- 
deed, from the experimental point of view, this system 
probably provides the best example of random quantum 
critical behavior. 

Theoretically, however, such systems still present a 
grand puzzle, in which but a few pieces are in place. Some 
analytical progress has been mack rfor metal-insulator 
transitions in d = 2 + e dimensions B E But in 2d, despite 
a set of simple non-interacting Fermion models which de- 
scribe the IQHE plateau transition, a controlled analytic 
treatment is sorely missing. Efficient numerical methods 
have been developed to investigate these models, andcro- 
vide a number of significant empirical observations JjI3 In 
both these cases, scaling is manifest in the vicinity of 
the critical point, with a diverging localization length 
£ ~ \M \~ v , where M measures the deviation from criti- 
cality. For the IQHE, v ss 7/3. For distances shorter than 



£, the single-particle electronic wavefunctions (i/K^W are 
extended, but exhibit complex multifractal scalinga. In 
particular, each disorder-averaged moment scales with an 
independent pair of critical exponents, which we denote 
as r(q) and y(q) for the q th moment: 
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Here the square brackets denote an ensemble average over 
disorder configurations, L is the linear extent of the sys- 
tem and \x\ <C L is assumed. The scaling for essentially 
all physical quantities can be formulated in terms of the 
set of exponents v,r(q), and y(q) (a simple example is 
given m Ref. 0). 

Supersymmetric (SUSY) techniques offer the tantaliz- 
ing possibility of a field-theoretic treatment of such de- 
localization transitions. SUSY has a long history in dis- 
ordered electronic systems, where it was first introduced 
by Efetov to describe diffusive metals E The SUSY non- 
linear sigma model, when linearized, provides a Gaussian 
or free field description of a diffusive metal. Expansions 
about the metallic phase in d = 2 + e dimensions give 
a fixed point which can be extrapolated to describe a 
3d metal-insulator transition. However, isolated der- 
ealization transitions in d < 2, such as in the IQHE, 
do not afford the luxury of expanding about a diffusive 
metallic phase. Recently, Zirnbauer and others have used 
SUSY to map the Chalker-Coddington network model for 
the IQHE transition into an interacting 1+1-dimensional 
field theory - a supersymmetric antiferromagnetic spin 
chain .El EiJ Unfortunately, this model has resisted all at- 
tempts at an analytic treatment, despite the essentially 
complete solution of a related supersymmetric ferromag- 
netic chaip-,which describes transport in a dirty 2d chi- 
ral metal.ta Falco and Efetov have recently applied the 
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SUSY non-linear sigma model in 2d to extract multifrac- 
tal wavefunction correlations, but in a crossover regime 
rather than at an isolated derealization transition.Ej 

Some analytic progress has been made using a toy 
model of Dirac fcrmions in a random vector potestta! 
(RVP), which exhibits a 2d derealization transition.liirtZl 
This model has the simplifying feature that a zero en- 
ergy (critical) wavefunction can be found exactly for 
any realization of the disorder, which enables analytic 
study of wavefunction multifractality. A number of dif- 
ferent formulations are possible, but a particularly in- 
triguing-approach involves mapping to 2d Liouville field 
theory,H3 which has been extensively studied in string 
theory. Away from criticality, however, far fewer results 
are known, and at present SUSY techniques have not 
been successful in this regard. 

In this paper, we study a one-dimensional tight- 
binding system of spinless fermions with random hopping 
matrix elements, 



H = - 



E 



^n(c n Cyj 
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Here the d s are canonical Fermion operators satisfying 
{cmjcjj} = S mn , and the random hopping strengths t n 
can be taken positive without loss of generality. The con- 
tinuum limit of this model is in fact a one-dimensional 
analog of the 2d RVP theory, and many of the same 
properties obtain. An exact zero energy wavefunction 
is known for each realization of the disorder. Critical 
singularities are present in the single-particle density of 
states as in the 2d RVP model (but in contrast to the 
IQHE transition). Indeed, for the Id random hopping 
model the density of states diverges at the band center. 

There has been considerable prior work on the Id tran- 
sition in the random hopping model, primarily focussing 
on properties derivable from the mean local Green's func- 
tion: the mean density of states and the typical local- 
ization length. This work was recently summarized in 
Ref. fill Employing a real-space RG method, D. S. Fisher 
(DSF) has obtained the spatial dependence of mean spin- 
correlation functions in several closely related Id models: 
the McCoy- Wu random transvere field Ising. model and 
random Heisenberg and XX spin chains .EjEj In this pa- 
per, we extend the above analyses using SUSY meth- 
ods to obtain the spatial dependence of the exact critical 
and off-critical scaling functions for t he mean F ermion 
Green's function, summarized in Eqs. ( 5.31 -5]4fl\— .(Un- 



published work by DSF using real-space methodscu cor- 
roborates our results). An important feature not present 
in the local properties is the emergence of a mean lo- 
calization length which controls the spatial decay of the 
average Green's function, 
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with e the energy from the center of the band. This 
length is much longer than the typical localization length 



£ e ~ | In e| , found previously by many authors. This im- 
portant distinction between typical and mean correlation 
lengths has been emphasized by D. S. Fisher in his anal- 
ysis of Id random spin-chains. 

In contrast to the IQHE transition which is mapped 
into a SUSY spin-chain, the SUSY formulation of the 
Id random hopping model is equivalent to the quantum 
mechanics of a single superspin, with "Hamiltonian" 



H = 2loJ z + 2m Q J x - 4g (J x y 
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The superspin operators J are defined in section III. This 
simplification enables us to systematically carry through 
the analysis from start to finish to obtain the exact crit- 
ical and off-critical scaling functions. A key motivation 
for doing this was to investigate in detail the novel el- 
ements which arise in a SUSY formulation of a der- 
ealization transition. Indeed, the "spin" Hamiltonian 
in Eq. (1.4) has many unconventional properties. It is 



non-Hcrmitian, requiring a distinction between left and 
right cigenstates. The spin operators themselves are in 
fact elements of the non-compact superalgebra U(2| 1,1). 
The non-compactness of the SU(1,1) subalgebra is man- 
ifest in the representations of the spin operators which 
are infinite-dimensional (i.e. the ladder of discrete J z 
eigenstates is infinite). As demonstrated in sections IV- 
V, near criticality the system explores the far reaches 
of this manifold of spin states, in a manner which can 
be described by Liouville quantum mechanics, which was 
recently introduced in studies of the zero-energy eigen- 
states in Ref. ^2|. Furthermore, the SUSY Hamiltonian 
is defective, i.e. the (right) eigenstates do not span the 
Hilbert space. To surmount this difficulty requires the 
definition of "pseudo-eigenstates" to complete the eigen- 
basis. A knowledge and familiarity of these features will 
likely be crucial to the success of future work applying 
SUSY to other (e.g. 2d) critical points. 

As a major function of this paper is pedagogy, we have 
attempted to present the material in enough detail to al- 
low the reader to appreciate the technical elements of the 
calculations. Section II discusses the model, its relations 
to various random spin chains, the continuum limit, and 
the relevant single-particle Green's functions. In section 
III, we describe the mapping to quantum mechanics, de- 
rive the SUSY Hamiltonian and its (super)symmetries, 
and detail the organization of states into superspin lad- 
ders and super-multiplets. The exact ground state and 
a class of excited states needed to compute the desired 
correlators are found in sections IV and V, leading to 
the final results in Eqs. (5.31 |5.4C ). Lastly, in section VI 
we pursue the Liouville quantum mechanics formulation, 
extending the treatment of Ref. ^2] to determine the full 
set of multifractal exponents 



r(q) = 0, y(q) = 3/2. 



(1.5) 
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II. THE MODEL AND SYMMETRIES 



A. Lattice model and continuum limit 



We begin with the free-fermion model, Eq. (|l.2| ). We 
assume that t n consists of a large uniform part, i, and a 
small random piece, St n . Under a Jordan Wigner trans- 
formation, this model is equivalent to a random exchange 
spin- 1/2 XX chain, 



Hxx — — ^ 2tn(SnSn + i + S„Sn +1 ), 



(2.1) 



where S n = <r n /2, with <x the usual vector of Pauli ma- 
trices. 

For uniform hopping the single particle states are 
plane waves, and H describes a band at half filling, 
with zero Fermi energy, and two Fermi points at kp± — 
±7r/2. With a small random component in the hopping 
strengths, the single particle states will be localized away 
from the band center, but due to a special particle hole 
symmetry (see below) the localization length diverges 
upon approaching zero energy. The density of states is 
also singular at zero energy. To study this derealization 
transition, it suffices to focus on states near zero energy. 
Provided St n is small compared to the band width t, it is 
legitimate to take a continuum limit, retaining a narrow 
shell of pure energy levels near the two Fermi points. 

To this end, we decompose the Fermi fields as 



(2.2) 



and assume that ip varies slowly with n. To take the con- 
tinuum limit we replace n by a continuous coordinate x, 
and approximate discrete differences with x-derivatives. 
For the Fermion hopping term this gives, 



c n c n+l 



-2i(-l)' 



(2.3) 



For uniform hopping, t n — t — > dx/2, the second term is 
rapidly varying and can be ignored, giving the expected 
(pure) Hamiltonian, 



Ho 



dx[tfj R id x t/j R 



iplidxi/jL,} 



(2.4) 



which describes right and left moving modes at the two 
Fermi points. 

A small random hopping 5t n causes scattering between 
the plane wave states. The important Fourier compo- 
nents of 8t n are at 7r, since these cause backscattering 
between the right and left movers. We thus decompose 
St n ~ * ( — l) n m(x)dx/2, where m(x) is assumed slowly 
varying. From the second term in (p. 3]), this leads to a 
(random) backscattering term in the continuum limit: 



Hi = -i f dx m{x) (^y L - ^Vfl) • (2.5) 



Employing a spinor notation, ip = (tJjr, iPl), the full con- 
tinuum Hamiltonian, H c = Ho + Hi takes the form 



H c 



dx ijy hip, 



with a single-particle Hamiltonian 

h = —i<j z d x + m(x)a y . 



(2.6) 



(2.7) 



It is convenient to decompose the function m(x) into 
a uniform and random piece as 



m(x) = mo + m(x), 



(2.8) 



where [m] cns . = 0, with the square brackets denoting an 
ensemble average. Non-zero mo correspondes to a (uni- 
form) staggering in the hopping, St n ~ (— l)"mo, which 
opens a gap in the pure spectrum about the band cen- 
ter. In the XX spin-chain, non-zero mo corresponds to a 
dimerization in the bond strengths, and the gap is a spin- 
gap due to singlet formation across the stronger bonds. 
Both mo and the energy e are tuning parameters which 
take one away from the delocalized critical point. 



B. Symmetries and derealization 



The lattice free Fermion Hamiltonian Eq. (1.2) is in- 
variant under the canonical transformation, 



(-1) 



(2.9) 



due to time-reversal and particle-hole symmetries, 
present even with random hopping strengths. As a con- 
sequence of this symmetry, the single-particle wave func- 
tions can be chosen real (time reversal invariance) and 
come in conjugate pairs with energy ±e. Specifically, for 
a given cigenfunction, 4> t (n) at energy e, there is a partner 
eigenstate with energy — e, given by 0_ e — (~l) n 4> e (n). 
At zero energy one thus anticipates special properties, as 
discussed below. 



In the continuum, the symmetry Eq. (2.E) becomes an 
invariance of H c under the canonical transformation, 



(2.10) 



for a = R, L. This symmetry restricts the allowed form of 
the single particle Hamiltonian, h. Specifically, h cannot 
have terms (with no gradients) proportional to a x , a z or 
the identity. A generic random Id tight binding Fermion 
model, in which the density of states is regular and all the 
eigenstates are localized, would not be particle-hole sym- 
metric, and additional terms, such as a spatially random 
a x term, would be present in the continuum Hamiltonian. 
The above symmetry is clearly crucial for the existence 
of derealization at the band center, e = 0. 
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An (un-normalized) extended state at zero energy can 
in fact be directly extracted from the continuum wave 
equation: h&(x) = 0, where $ is a two-component wave 
function. Writing 



$(x) = <j>±{x) 



1 

±1 



the scalar function <p(x) satisifies, 

(d x ± m(x))(j>± = 0. 
This can be integrated to give 

/ / \ ±1 dx'm(x') 

<p±{X) (X e J '. 



(2.11) 



(2.12) 



(2.13) 



For random m(x), with mean zero, this wavefunction 
is clearly not exponentially localized. If the random 
function m(x) has short-ranged spatial correlations, the 
logarithm of the wave function undergoes a Id random 
walk. For a Gaussian distribution of m(x) the (unnormal- 
ized) wavefunction is log-normally distributed. This wave 
function is a one-dimensional analog of the exact zero 
energy wavefunctions written down for 2d free Fermions 
described by a Dirac theory with random vector poten- 
tial. As in the 2d case, the wavefunction is very broadly 
distributed, and it's correlations can be characterized by 
a multi-fractal scaling description. We return to a discus- 
sion of the multi-fractal characteristics of this wave func- 



tion in Section VI, where we compute the multi- fractal 
spectruipiexplicitly, following recent work by Shelton and 
Tsvelik.E3 

Away from criticality, f or no n-zero mo, the zero en- 
ergy wavefunctions in Eq. ( [2.13] ) are exponentially grow- 
ing and decaying functions, (f>±(x) ~ e ±m o x _ While they 
are non-normalizablc in infinite space, for a finite system 
they describe solutions which decay exponentially into 
the system, with an associated localization (or correla- 
tion) length, £ = l/mg. The critical exponent, defined 



(2.14) 



is v = 1. As emphasized by DSF, in addition to the 
length £ which describes the decay of a typical (unaver- 
aged) correlation function, there is another divergent 
length, £, which describes the decay of ensemble aver- 
aged correlation functions. Consistent with arguments 
by DSF, we find below that this latter length diverges 
more rapidly with an exponent v = 2. 

Two similar lengths may be defined by approaching 
the critical point at finite energy e, but with zero mass 
mo = 0. Using the Thouless construction from the local 
Green's function, previous authors have found a (typical) 
localization length £ c ~ |lne|. In constrast, employing 
a real-space RG approach, DSF has shown that mean 
correlation functions decay with a longer length, which 



Another important characteristic of the above exact 
zero energy wavefunction is that it is nodeless, for each 
and every realization of the random potential m(x). Be- 
cause of this, critical properties of the Id localization 
transition at e = are contained in the ensemble aver- 
aged single-particle Green's function, in contrast to the 
conventional Anderson transition in higher dimensions. 

Below we briefly consider symmetry properties of the 
single Fcrmion Green's function, and obtain expressions 
in the continuum limit. The next sections are devoted to 
evaluating the ensemble averaged Green's function using 
supersymmetry methods. 



C. Green's Functions 



Consider the single Fermion Green's function at energy 
e, defined as 



G(n,n';e + iLu)=i / dte^+^^Mi)^, (0)|t>), 



(2.15) 

where \v) denotes the Fermion vacuum, and c(t) = 
e lUt ce~ lHt with Ti the lattice Hamiltonian. Here u> is 
a small imaginary part to the energy. In practice below, 
we will take the real part of the energy to be zero, e = 0, 
calculate Q{iui) for real ui, and then extract the energy 
dependence via an analytic continuation. 

The spectral decomposition of Q in terms of the exact 
eigenstates, <j> e , takes the form 



g(n,n';iw) = ^ 



4> e {n)4>e{n') 



(2.16) 



Using the symmetry property 0_ e (n) = (— l)™0 £ (n), one 
can readily show that Q(n,n' ,ico) is real and even in uj 
for n — n' odd, and purely imaginary and odd in u fo r 
n — n' even. For example, with n — n' even, Eq. (2T£) 
can be rewritten as 



Q(n, n'; iuj) 



£>0 



(2.17) 



The Green's function, as defined in Eq. (2.15), can be 
re-expressed in terms of the continuum Fermion fields, 



ijj(x), by employing the decomposition Eq. (2.2). For 



varies as £ e 



In ( 



n — n' 3> 1, the discrete separation can be replaced by a 
continuous distance x — x' (lattice constant equal to one). 
One finds 

g(n, n'; iuj) = ^(-l^+^'G^, x'; iu), (2.18) 

a, 13 

with a and (3 running over the two spinor components, 
denoted as either (R,L), (0, 1) or (t,|)- Here G a p is de- 
fined in terms of the continuum Fcrmion fields, ip, as 
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G a p(x, x'; iuS) = i 



dte- ut (v\i/j a (x,t)ijl(x' ,0)|»), 

(2.19) 



with ip(t) = e lWc *^/;e~ lWct and 7i c the continuum Hamil- 
tonian. This continuum Green's function can alterna- 
tively be express ed in terms of the single-particle Hamil- 
tonian, h in Eq. (2/7), as 



G a p(x,x';iu>) — (x,a\- 



1 



\x',p), 



(2.20) 



where \x, a) denotes a Fermion at position x with "spin"- 
component a. 

Of interest is the behavior of the single-particle Green's 
function upon ensemble averaging over disorder realiza- 
tions. To be concrete, we take the random function m(x) 
to be Gaussian with, 



[rh(x)rh(x')] ens . = 2gS(x — x'). 



(2.21) 



Ensemble averaged Green's functions, which we denote 
with an overbar, become translationally invariant: 



1 N 

G{x;iu) = — 22[G(n + x,n;iuj)] ens ., 



(2.22) 



with N — > oo the numb er of sites in the tight binding 
lattice. From Eq. ( 2.18 ), this is related to the averaged 
continuum Green's functions, 

G a/3 (x; ioj) = [G a p{x, 0; w)] ens ., (2.23) 



via 



G(x; iu) = i x 2j(-l) Q "G QQ (x; iw), (2.24) 

OL 

with the separation x either even or odd. 

The mean density of states for the original lattice 
Fermions can be written, 



p(e) = lim — ImQ(x = 0; e + iu>). 



(2.25) 



We shall also be interested in the spatial dependence of 
the correlation function, 



C(x,e) = -lmG(x;e + iO + ). 

TT 



(2.26) 



This function is expected to decay exponentially with a 
mean correlation length £. Due to the delocalized zero 
energy wave function (at mo = 0), £(e) should diverge 
upon approaching the band center, e — > 0. 

In Section III we will construct a generating functional 
which can be used to extract G. Our strategy will be to 
calculate G(x,iu>) for real u>, and then perform an an- 
alytic continuation to extract the density of states and 
C(x,e). 



D. Related Random models 



As noted above, under a Jordan- Wigner transforma- 
tion the lattice free- Fermion Hamiltonian H in Eq. (1.2) 
is identical to a random exchange spin-1/2 XX chain. 
Some properties of the spin-chain can be extracted from 
the Fermion density of states, specifically the specific 
heat and the z-component magnetization, (S z ), in re- 
sponse to a magnetic field along the z-axis. Unfortu- 
nately, spin correlation functions are notoriously difficult 
to extract from the Free-fermion representation, due to 
the non-local relation between spin and Fermion opera- 
tors (Jordan- Wigner string). Nevertheless, one expects 
that the correlation decay length of the Fermion Green's 
functions will also control the decay of spin correlations. 

As shown originally by Shankar and Murthy, a second 
spin model which is equivalent upon fermionization to 
the free Fermion model TL, is the Id random quantum 
Ising chain in transverse field, with Hamiltonian 



n 1 = 



(2.27) 



Here K\ and Ki are spatially random field and Ising 
exchange constants, respectively. This model is an 
anisotropic ( "time-continuum" ) version of a 2d classical 
Ising model with random exchange interactions perfectly 
correlated in one of the two directions, a model first stud- 
ied by McCoy and Wu. For completeness, we sketch this 
fermionization procedure in Appendix A, where we show 
that the low energy properties of Hi follow from the prop- 
erties of the single-particle Hamiltonian h in Eq. (2.7). 

Finally, we should mention the equivalence between the 
free Fermion Hamiltonian H. and a Id model of quan- 
tum particles connected by random strength harmonic 
springs, a model firstintroduced and analyzed by Dyson 
almost 50 years ago.E3 



III. QUANTUM MECHANICS 



As discussed above, much of the information of interest 
is contained in the mean single-particle Green's function 
G a p{x, e + iu). In the following we will construct a gen- 
erating functional which can be used to extract G. Our 
strategy will be to calculate G(x, iuj) for realuj, and then 
perform an analytic continuation to extract the density 
of states and C(x,e). 
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A. SUSY generating functional 



lim STre- iH , 



(3.6) 



Our analysis is based on employing the well-known 
field-theoretic representation of an operator inverse to 
express the ensemble averaged Green's function, 



G a p(x,iw) = [< x,a\- 



1 



■\O,0>] 



(3.1) 



(3.2) 



(3.3) 



Here the functional integration is over Grassmann fields 
ij){x) and complex fields £(x), with the action 



h — iuj 
as a functional integral, 

G a p(x,iuj) = i(ip a (x)ipp(0))s. 

with 

(0) s = [ DiPD~$DiDCOe- s 



S= dx$(ih + w)i> + £*(ih + tj)ti]. 



(3.4) 



The most noteworthy point is that in Eq. ( |3.3|) , we have 
included the complex scalar field £ in order to cancel the 
Fermionic determinant which naturally occurs due to the 
if) integratio n. In doing so, we obtain the supersymmetric 
form in Eq. (3.4). This enables an ensemble average over 
the Gaussian disorder rh(x) to be readily performed. 

For ease in presentation, it is convenient to speak of a 
slightly simpler object, the partition function 



Z = / DipDipD£D£* 



(3.5) 



where the functional integration is, as before, over Grass- 
mann field s ip(x) and complex fields £(x), with the action 
of Eq. (3.4). Correlation functions are obtained by simply 
inserting the appropriate fields after the integration mea- 
sure and ensemble averaging. The crucial cancellation of 
the Fermionic and Bosonic determinants then gives the 
trivial identity (1) = 1, or Z = 1. From the functional 
integral formulation, the reason for generating Green's 
functions for real oj becomes clear: while the Fermionic 
functional integral is always well defined, the bosonic one 
is only convergent provided that the action is bounded 
below. This is the case here provided u> > 0, since h is 
hcrmitian. 



B. Tranformation to quantum formulation 

The above action corresponds to a random one- 
dimensional statistical mechanics problem. After ensem- 
ble averaging over the random function m(x), the model 
is translationally invariant. Our approach is to extract 
the transfer matrix, T = e~ H , which can be used to re- 
construct the averaged generating function: 



where L is the length of the system. The symbol STr 
indicates the supertrace, defined by 



STrO = Tr [(-l) N fQ] , 



(3.7) 



where Tr is the conventional trace, and Nf is a Fermion 
number operator defined below. Although H is an opera- 
tor, easily expressed in terms of Fermi and Bose operators 
(see below), it is "zero-dimensional," being independent 
of the spatial coordinate x. The problem is thus reduced 
to studying the quantum mechanics of the "Hamiltonian" 
H. As usual, the spectrum of H contains information 
about the correlation length of the Id system - in this 
case the localization length. 

To extract the operator H, we massage the action 
Eq. (3.4) into the form of a coherent state path integral 
with x playing the role of imaginary time. To this end we 
let £| — > — £|> leaving £j unchanged. Similarly we trans- 
form the independent Grassmann fields as ip± — > — -0^ 
and 0i — > ipi, leaving the spin-up Grassmann fields un- 
changed. The action can then be written 



S = dx(C +£>u + C m ), 



with 



C = ipd x ip + C^xi, 



and C r , 



m{x)A with 



ip+ip, - ip^ipi + (ip ->■ f). 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



Notice that Co is now in the standard form for a coherent 
state path integral if x is reinterpreted as an imaginary 
time coordinate. 

Before extracting H, we perform an ensemble average 
over the disorder. Since we have assumed a Gaussian 
distribution, this is readily performed to extract [i?] ons .. 
The only term in the action which is modified is C m , 
which now becomes, 



CZ S - =moA-gA 2 



(3.12) 



The transfer "Hamiltonian" H can now be read off, 
since the full action takes the form, S = f (Co+H(ip, £)). 
In passing to the Hamiltonian framework, the Grassmann 
fields are replaced by Fermion operators, tp — > /, ip — > f\ 
and the complex fields by Bose operators, ^ — > £>, ^* — > b\ 
where / and b satisfy canonical commutation relations: 



[fa, fp]- = [b a ,b\] = 5 afi . 
The resulting "Hamiltonian" is, 



(3.13) 
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H = cu[pf + 6 f 6] + m A - gA 2 , 



(3.14) 



C. Supersymmetry 



with 



.4 



fUl ~ frh + (/ - b). 



(3.15) 



Although we will hereafter refer to H as a Hamiltonian, 
it is important to keep in mind that this operator is not 
Hermitian. 

Since H does not conserve the Fermion number, Nf = 
f'f, it is convenient to perform a particle hole trans- 
formation, defining new Fermion fields via a canonical 
transformation 

*r = /r. n = /|. (3-16) 

where F^F commutes with H . To preserve the Bosc- 
Fermi supersymmetry one can also define, 



B i 



Si 



(3.17) 



However, i?| does not satisfy the canonical Boson com- 
mutator, but rather, = — 1. To restore the 
canonical form we define 



so that B and B satisfy, 

[B a ,Bp] = 5 a p. 



(3.18) 



(3.19) 



However, it must be kept in mind that B ^ B^ . 

In term of these new operators, the Hamiltonian H 
becomes 



H = Lo[F^a z F + Bcj z B]+moA-gA 2 , (3.20) 



with 



A = F^a x F + Ba x B. (3.21) 



At this sta ge it is convenient to express the Green's 
function Eq. ( ]3 . 1[ ) as a supertrace ov er q uantum states. 
Consider the representation Eq. (3.2) in terms of 
Fermions. After transforming to the coherent state path 
integral form, the Grassmann fields can be replaced by 
Fer mion operators: ip(x) — ► e fe~ xH . The average in 
Eq. (3.2) is replaced by a supertrace over quantum states: 

(Oty,0) S -> (O(f.b)) = STr[Oe- LH }. (3.22) 

This gives, 

ZS a p(x;iu) = z(-l) a STr[i^(:r)i^ e - iH ], (3.23) 
with F a (x) = e xH F a e- xH . 



In order to disucss the symmetries of the Hamiltonian 
H it is useful to introduce a four-component superfield, 



(3.24) 



We will use latin indices (a, &,...) to denote the 
Fermion/Boson label, i.e. a = B,F 0, 1. From this 
superfield, one may build a three-component superspin, 



3 'ab= 2*a" Cr a/3^'b/3, 



(3.25) 



where sums on the greek indices are implied. One can 
also define a set of "charges" , 



Qab = ^aa^ba, 



(3.26) 



(sum on a). 

The diagonal components of 3 have special signifi- 
cance. In the Fermionic sector, 



J ll = S = F^-F 



(3.27) 



forms a set of ordinary Hermitian SU(2) spin operators, 
satisfying [S'%<S J ] = ie 13 S . In the bosonic sector, we 
may similarly define three other currents, 



3 



oo 



J = B°-B. 



(3.28) 



These also satisfy the SU(2) algebra, [J\P\ = ie ijk J k . 
However, although J z is Hermitian, J x and J v are antir 
Hermitian. One could define a Hermitian set of opera- 
tors, multiplying the x and y components of J by a factor 
of i. These would then satisfy SU(1,1) commutation re- 
lations, instead of SU(2). 

A useful object is the total spin-current 



3 = 3 a 



(3.29) 



where again the sum on repeated indices is implied. This 
spin-current commutes with the charges: 



[J, Qab] = 0. 



(3.30) 



Since the Hamiltonian can be expressed in terms of this 
total spin-current, 



H = 2uJ z + 2m J x - Ag (J x f , 



(3.31) 



the charges Q ao also commute with H. Thus Q ao gen- 
erate a set of (super)symmetries of H . Because of their 
importance, it is convenient to name them individually: 

N B = Q 00 = BB, N F = Q n = (3.32) 
Q = Qoi=BF, Q = Q W =F^B. (3.33) 

These latter two operators are Fermionic "charges", 
which will be extremely useful in determining the ground 
state of H. They obey 
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Q 2 = Q = 0, 



(3.34) 
(3.35) 



{Q,Q} = N, 

where the total charge is defined as 

N = N B + N F , (3.36) 
which commutes with all sixteen of the U(2|l,l) currents. 

D. Hilbert space and Representations 

Finding the ground state and low energy excitations 
of H is complicated by the enormity of the Hilbert 
space. Indeed, since the number of bosons with spin 
a, n a — b' a b a , is unbounded, the Hilbert space is ac- 
tually infinite. One basis of states spanning the Hilbert 
space may be written as a direct product of bosonic and 
Fcrmionic states: 



|n-|-nj.a) = |n| n i) ® 



(3.37) 



where n|,ri| = O...00 are the number of up- and down- 
spin bosons, respectively. The parameter a F labels the 
Fcrmionic sector, which is spanned by the Fcrmionic vac- 
uum I vac), which is anihillated by F, and three other 
states 



I TH-F-flvac), 
I |) = F/|vac), 
|T> = i^/|vac>. 



(3.38) 
(3.39) 
(3.40) 



As is usual in a quantum mechanics problem, we may 
simplify matters by choosing a maximal set of commut- 
ing variables, whose eigenvalues are good quantum num- 
bers. In this case, Nb and N F are obvious choices. The 
Fermionic charges Q and Q cannot, of course, be diag- 
onalized. They can, however, be combined to form two 
additional bosonic charges 



r = QQ, r = QQ. 
It is straightforward to show that 



T 2 = AT, 

r + r = a, 



r = 
rr 



NT, 
-- 0. 



(3.41) 



(3.42) 
(3.43) 



These relations imply that the eigenvalues of (r, r) are 
either (0, A) or (A,0). The four operators N B ,N F ,T,T 
form the desired set of good quantum numbers, and it is 
desirable to reorganize the states given above into a basis 
diagonal in these variables. 



1. SU(2) 

Before proceeding to determine this basis, consider 
first the fermionic sector of the theory. The fermion num- 
ber is in fact related to the total spin via 



with 



S 2 =8(8 + 1), 



s = N F (2-N F )/2. 



(3.44) 



(3.45) 



So we can think of N F as determining the representation 
of SU(2). Note that the singlet (s = 0) representation 
occurs twice - for N F = 0, 2. 



2. SU(1,1) and bosonic ladders 

Similarly, the bosonic states may be separated into 
multiplets with fixed Nb = rif — — 1. Each such mul- 
tiplet is actually a distinct irreducible representation of 
SU(1,1). To see this, consider the Casimir operator 

j2 = (jx)2 + (jy)2 + (J , )2 = ( N 2 + 2N B )/4. (3.46) 

Fixing N B thus fixes the "total spin" of the SU(1,1) rep- 
resentation. Following the analogy with SU(2), we may 
label the states by their total spin and, e.g. the spin J z 
along the z-axis, 



J 2 \jn)=j(j + l)\jn), 
N B \jn) = 2j\jn), 

j-w- ri + w+11 



+ n 



\jn). 



(3.47) 
(3.48) 

(3.49) 



Note that in the last equation we have departed from 
the usual convention for denoting J z eigenvalues. This 
is convenient because the quantum number n as we have 
defined it takes integer values n = 0, . . .00. The total 
spin can take half integer values j = 0, ±1/2, ±1, . . .. 

It is also helpful to have explicit expressions in terms 
of the previous basis. Since the j th block of states cor- 
responds to a ladder satisfying rif — = 2j + 1, it can 
be conveniently generated using the raising and lowering 
operators: J ± = J x ± iJ v . As usual, the lowest weight 
state in each ladder, denoted \j0), is constructed to be 
annihilated by J_ = — 



2j+l 



(3.50) 



where \v) denotes the bosonic vacuum: b a \v) — 0. Each 
ladder is constructed by acting on the associated lowest 
weight state with powers of J+ = fejbj; 



j + l\\ 



v n!(n +|2j + l|)! 
where n runs from zero to 00. 



1/2 



(J + ) n \j0), (3.51) 
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3. SUSY ladders 

Clearly the set of eigenvalues of Nb, Np, T, T is insuf- 
ficient to distinguish all the states in the Hilbert space. 
To provide a unique labelling, we will choose to diago- 
nalize the additional operator J z = J z + S z . This choice 
is natural in that J z commutes with the other four di- 
agonal charges. Moreover, the other current present in 
the Hamiltonian, J x , leaves the original four quantum 
numbers unchanged, mixing only different values of J z . 
The collection of states with different J z but with the 
other four charges fixed may be viewed as the basis for a 



representation of the algebra of the J7" operators, i.e. a 
peculiar (non-unitary) representation of SU(2). 

Two sets of such states are easily constructed. These 
are 

|2j, 0, 0, N, (\N + 1| + l)/2 + n) = \jn) ® |vac), (3.52) 

\2j, 2, N, 0, (|JV - 1| + l)/2 + n) = \jn) ® | |T>, (3.53) 

where the quantum numbers inside the bras on the left- 
hand-side denote eigenvalues of (N B , Np, T, T, J z ), in 
that order. 



States with Np = 1 are slightly more complicated, since they can involve linear combinations of up or down 
fermions. For N ^ 0, these are 



|2j, 1, (N - |AT|)/2, (N + |7V|)/2, \N\/2 + n) = J £±W\ jn ) ® | |) + J^JL-\j n - 1) ® | T >, (3.54) 



\2j, 1, (N + \N\)/2, (N - \N\)/2, \N\/2 + n) = ^^-2_|j„) ® | |) + - 1) ® | T), (3.55) 

where in the first set above n = 0, 1, 2, . . . oo, while in the second set n = 1, 2, . . . oo. 

Apparently these two ladders of states become identical for N = 0. What happens in that case? There are two 
states which are eigenstates of Nb = —l,Np = 1, and J z = n: 

| - 1/2, n) ® | |), | - 1/2, n- 1)® |T). (3.56) 

If, however, we attempt to diagonalize, e.g. T, in this two-dimensional space, we find that it is impossible! In fact, 
there is only a single eigenstate for each n, 

|_ M)0 , ,n) = (^[|-lA«>®U> + |-lAn-l>®IT>] n>0 

1 ' \ | - 1/2,0) ® | |) n = V ; 

which are annihilated by both Q and Q. Clearly this set of states does not span the full subspacc. To complete it, 
we may define an orthogonal ladder of states 

|-1,1, *, *, n) = ±= [| - 1/2, n) ® | |) - | - 1/2, n - 1) ® | |)] , (3.58) 

for n = 1,2,... oo. It is important to note that | — 1,1, *, *, n) is not an eigenstate ofT and T. Instead, acting with 
these operators on | — 1, 1, *, *, n) gives back | — 1,1, 0, 0, n), i.e. T and T act like projection operators in this subspace. 



E. Eigenstates 

1. supermultiplets 

In a conventional quantum mechanics problem, we can 
look for eigenstates of the Hamiltonian separately for 
each distinct set of eigenvalues of the chosen commut- 
ing operators. As we have seen above, most of the states 
in the Hilbert space can be specified in this way. To avoid 
the exceptions for the moment, consider first the sectors 
with N 7^ 0. Then we may find right eigenstates of the 
Hamiltonian 

H\N B ,N F ,r,T,E) =E\N B ,Np,T,T,E), (3.59) 



where each such state can be expanded in a basis of ap- 
propriate J z eigenstates, i.e. 

\n b ,n f , r, r, E)=J2 xe b ' Nf ^ y (J z )\Nb,Np,t,T, j z ). 

J* 

(3.60) 

In a theory with ordinary bosonic symmetries, com- 
pletely specifying the quantum numbers usually deter- 
mines a unique set of energies - accidental degeneracies 
are rare. In a SUSY theory, however, the additional 
fcrmionic generators Q and Q lead to extra relations be- 
tween states. 
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To see this, note that since Q and Q commute with H, 
acting with them upon an eigenstate must either produce 
another eigenstate with the same energy or annihilate it. 
Because Q and Q do not commute with the four diago- 
nal operators, however, they must change these eigenval- 
ues, and hence connect distinct states. Because the total 
charge N does commute with Q and Q, it (along with 
the energy E) can be used to characterize a multiplet of 
states connected in this way. 



l,l,0,0,n) 



(3.68) 



N, 



1 



N = 



10) 





(b) 



Fig. 1: Organization of multiplets in the SUSY Hamil- 
tonian. (a): Doublets, involving one bosonic and one 
Fermionic state for total number N ^ 0. (b) Singlet 
ground state |0) and quadruplets, with two Fermionic 
and two bosonic states for N = 0. See text for details. 



and is hence orthogonal to | — 1, 1,*, *,n). 

Taken together with the other three states of energy 
E, | — 1, 1, *, *, E) forms the fourth mem ber of a super- 
quadruplet (see Fig. lb). Using Eq. ( 3.67 ), it is straight- 
forward to see that the three true eigenstates may be ob- 
tained through the action of Q and Q on | — 1, 1, *, *, E), 
In particular, 

Q\ - 1, 1, *, *, E) = ci |0, 0, 0, 0, E), (3.69) 

Q\ - 1, 1, *, *,E) = c 2 \- 2, 2, 0, 0, E), (3.70) 

Q\0,0,0,0,E) = c 3 \- 1,1,0,0, E), (3.71) 

Q\ - 2, 2, 0, 0,E) = c A \- 1, 1, 0, 0, E), (3.72) 
Q\ -2, 2,0,0, E) = Q\0, 0,0,0, E) 

= Q\-1, 1, 0, 0, E) = Q\- 1, 1, 0, 0, E) = 0. (3.73) 

We expect the exception to this rule to be a unique 
ground state with zero energy, forming a super-singlet. 
We can essentially pin down the sector in which such 
a singlet can occur. It cannot occur for N ^ 0, since 
r + r = N ^ 0, and hence either Q or Q would create a 
new degenerate state. For N = 0, it is possible only for 
a state in the N F = -N B = 1, T = T = ladder. This 
is consistent with the existence of one additional state 
| — 1, 1, 0, 0, n = 0) in this basis. 



For N ^ 0, the states are in fact organized into dou- 
blets, as indicated graphically in Fig. la. One class of 
doublets includes the ladder of states with Np = 0, 

Q\N, 0, 0, N, E) = c\N -1,1, N,0,E), (3.61) 
Q\N-1,1,N,0,E) = {N/c)\N,0,0,N,E), (3.62) 
Q\N, 0, 0, N, E) = Q\N - 1, 1, N, 0, E) = 0, (3.63) 

while another includes the states with Np = 2, 

Q\N- 2,2,N,0,E) = c\N- 1,1, 0,N, E), (3.64) 
Q\N - 1,1,0, N,E) = (N/c)\N -2,2,N,0,E), (3.65) 
Q\N - 2, 2, N, 0, E) = Q\N - 1, 1, 0, N, E) = 0, (3.66) 

where different constants c can appear in different ex- 
pressions. 

For N = 0, things are more complicated. Distinct 
eigenstates can be found in the sectors with Np = 
~Nb — 0, 2, and for the ladder of states with Nf — 
—Nb = 1 and T = T = 0. The fourth subspace, spanned 
by the | — 1, 1, *, *, n) basis, however, is not closed under 
the action of H. The best that one can hope to achieve 
in this sector, therefore, is to find an eigenstate of the 
Hamiltonian projected back onto the same sector. Such 
a state | — 1,1,*,*, E) satisfies 

H\ - 1, 1, *, *, E) = E\- 1, 1, *, *, E) + (3.67) 

where 



2. partition function and non-hermiticity 

Thus far, we have described the organization of eigen- 
states. To calculate physical quantities, we will need to 
perform a trace over the Hilbert space of e~ LH . Once the 
Hamiltonian is diagonalized, this trace can be performed 
separately in each multiplet. Some care must be taken 
to account, however, for the non-Hermiticity of H. 

To do so, let us consider the behavior of the currents 
under hermitian conjugation. The usual SU(2) spin gen- 
erators S are of course Hermitian. The bosonic currents 
J are, however, mixed. In particular J z is hermitian 
while J x and J v are anti-Hermitian. Since the Hamilto- 
nian involves only J z and J x , the net effect of conju- 
gation is to change the sign of J x . This can in fact be 
accomplished by the unitary transformation (U = U') 

U = e™ r , (3.74) 

which rotates J x 

UU X U = -J x (3.75) 
and therefore conjugates the Hamiltonian 

U^HU = (3.76) 
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It therefore follows that if \E) R is a right eigenstate of H H^\E) L = E\E) L & L (E\ H = L (E\ E. (3.78) 

with energy E, 

Within each of the doublets, it is advantageous to nor- 
(3.77) malize the states such that 



\E) L =cU\E) R 
is a left-eigenstate of H, i.e. 



AE\ E)r = 1. 



(3.79) 



Left and right eigenstates with different quantum numbers are automatically orthogonal, so this normalization 
guarantees the resolutions of the identity 



IN: 



N:N F =0,l 



IN: 



N:N F = 1.2 



E 

E 

E 



\N, 0, 0, N, E) R L {N, 0, 0, N, E\+\N- 1, 1, N, 0, E) R l{N-1, 1, N, 0, E\ 
\N - 2, 2, N, 0, E) R L {N-2, 2, N, 0, E\+\N- 1, 1, 0, N, E) R L {N-1, 1,0, N, E\ 
for N 0. Inserting these into the usual expression for the trace in terms of an orthonormal basis gives 
STt n . Nf=0A O = J2 ~ l{N,0,Q,N,E\ O\N,0, 0,N,E) R + L {N-l, 1,N,0,E\ 0\N - 1, 1, N, 0, E) R 

E - 

STr A r ;A r F=1 , 2 (9 = J2 - l(N - 2, 2, TV, 0, E\ 0\N - 2, 2, N, 0, E) R + L (N - 1, 1, 0, TV, £| (9|7V - 1, 1, 0, N, E) R 



(3.80) 
(3.81) 

(3.82) 
, (3.83) 



where the minus signs arise from the (—1)^ in the def- 
inition of the supertrace, and again this holds only for 
N 7^ 0. Similarly, the ground state gives a unique contri- 
bution 



STr O = £<0| O\0) R , 



(3.84) 



where |0) = | — 1, 1, 0,0, E = 0) is the ground state. 

The quadruplets are somewhat trickier. The two states 
with s — (i.e. Np = 0, 2) in each quadruplet can be 
separately orthogonalized, i.e. 

L (0,0,0,0,£| 0,0,0,0, E) r = 1, (3.85) 
L {-2, 2, 0,0,E\ -2,2,0,0,E) R = 1, (3.86) 

since U\ — 2, 2, 0, 0, E) R remains in the subspace with 
N F = —N R = 2, and similarly for the Np = N R = 

The corresponding resolution of the identity is 



state. The same is not true for the two states with 
Np = —N R = 1- In fact, it is straightforward to show 
that corresponding left and right eigenstates are actually 
orthogonal! That is, 

L (-1,1,0,0,E| -1,1, 0,0,^ = 0, (3.87) 
L {-1,1,*,*,E\ -1,1,*,*,E) R = 0, (3.88) 

for E 7^ 0. This is because the U operator takes each 
basis state from one sub space into the other, as can be 
seen directly from Eqs. (3.57,3.58). Instead of the usual 
normalization condition, therefore, we must require 



L (-1,1,0,0,E\ -1,1,*,*,E) R =1, 
L (-1,1,*,*,E\ -1,1,0,0,^)^ = 1. 



(3.89) 
(3.90) 



1jv=o = 



\2,-2,0,0,E) RL (2,0,0,0,E\ + |0, 0, 0, 0, E) R L (0, 0, 0, 0, E\ 
1,1, 0,0,E) RL {-1,1, *,*,E\ +| -1,1, *,*, E) RL (-1, 1,0,0, E\ 



(3.91) 



From this identity follows the expression for the trace, 



STrjv^oO = 



L (2, -2, 0,0,^| O\2,-2,0,0,E) R - L (0,0,0,0,E\ 010,0,0,0,^ 



L {-1, 1,0, 0,E\ O|-l,l,*,*,^B+i(-l,l,*,*,f7|O|-l,l,0,0,£)j 



(3.92) 
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In this paper we are studying primarily properties of 
the system in the thermodynamic limit L — > oo. In this 
case the (super)trace of e~ LH reduces to an expectation 
value in the singlet ground state, i.e. 

STrOe _iH 



of SU(2) algebra of the total current As such, it must 
be annihilated by both Q and Q, and hence fits into the 
subspace with Np — —Nr = 1, T = T = 0, i.e. 



(O) 



,(0| O\0) R . 



(3.93) 



|0)k = X)^|- 1,1,0,0/, 



(4.1) 



It is worth noting, however, that SUSY indeed implies 
that the partition function itself is exactly one, even in 
the finite system. This is because each multiplet other 
than the singlet contains equal numbers of fermionic (Nf 
odd) and bosonic (Nf even) states, which thereby can- 
cel in the supertrace. All that remains is the expectation 
value in the zero-energy ground state, which is unity even 
for L finite. 



IV. GROUND STATE PROPERTIES 

A. Schrodinger Equation 

As we saw in the previous section, the ground state 
resides in a unique (non-unitary) singlet representation 



n=0 



In this basis, J x and J z have the very simple matrix 
elements 



J z \n} = n\n} , 
1 



(n + l)|n + l) -(n-l)|n-l)o 



(4.2) 
, (4.3) 



where we have introduced the more compact notation 
\n)o = | — 1,1,0,0, n). Using these results, and the re- 
quirement of a zero energy ground state, 

H\0) R = E\0) R = 0, (4.4) 

we obtain straightforwardly the Schrodinger equation 



(n + l)(p n +2 + 2ri(f> n - (n - l)4> n -2 ~ m (4> 



71+1 



4>n-l) + 2w</>„ 



= 0. 



Which implies 

- (n + l)(/) n+ 2 + 2n<p n - (n - l)(f> n -2 - 2M ((f>„+i - <j> n -i) + ^4>n = 0, for n ^ 0. 
Here we have defined M — mo/2 to simplify some expressions in what follows. 

B. Normalization and density of states 



(4.5) 



(4.6) 



where the second equality follows from Eq. (3.23). Simple 
manipulations then give 



From the previous section, given the right eigenstate 
|0)i?, a corresponding left eigenstate is obtained by 



5(0;^ W )=^£(-1^| 



(4.11) 



\Q)l=cU\Q) r = {-1) jz - 1 / 2 \Q) r , 



n=0 



(4.7) 



with an appropriate choice for the constant c. Applying 
this to the basis of J z eigenstates gives 



(-l)' r - 1/2 |->o = { 



_/(-l) n | -1,1,*, *,n) n>0 



n = 



(4.8) 



C. Continuum limit and boundary conditions 

In order to solve for the ground state wavefunction, 
first note the following simple property: if (j)^(M) is a 
solution of Eq. (4.6), then so is 



Since the \n} and | — 1,1, *, *, n) states are orthogonal, 
the normalization requirement reduces to 



<f>-(M) = (-ir<j>+(-M). 



(4.12) 



0)r = 



= 1. 



(4.9) 



T he de nsity of states is obtained in a similar way. From 
Eq. (2.25) the density of states is obtained by analytically 
continuing, 



G{x = 0; iuj) = G aa (0, 



Uxl 



-2i 



S Z \0) R , (4.10) 



Since this is a fourth-order linear difference equation, we 
would expect four linearly independent solutions. The 
above result reduces these to two trivially related pairs. 
To understand the nature of these two solutions, it is in- 
structive to consider the limit M = to = 0. In this case 
there are two obvious solutions: <f> n — 1 and (f> n = (— l) n 
satisfy the E = Schrodinger equation ( |4.6| ). However, 
these wavefunction are not normalizable, a feature due 
to a pathology of the theory at lu = 0. Indeed, non-zero 
uj is essential in guaranteeing convergence of the bosonic 
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functional integral in the generating function. Moreover, 
we are interested in behavior off criticality with M =/= 
and at finite energy iui — > e. 

A ge neral solution of the Schrodinger equation 
Eq. (4.6) is daunting (although possible for u> ^ - see 
Appendix B), so for simplicity we focus our attention on 
the critical regime very close to the localization transition 
where both u>, M <C g = 1. In this scaling regime one ex- 
pects the wavefunction to remain close to a superposition 
of the two trivial (constant and (—1)™) solutions, i.e. 



4> n = c 1( f>(n, M) + c 2 (-l) n 0(n,-M), 



(4.13) 



where cf>(n,M) is slowly varying with \<f>(n + 1,M) — 
4>(n,M)\ <C \<f)(n,M)\. This suggests a "continuum" 
limit, in which </>(n, M) may be regarded as a continuous 
function of n, and discrete differences in the Schrodinger 
equation are replaced with derivatives. In this contin- 
uum approximation the Schrodinger equation Eq. (4.6) 
becomes, 



+ (1 + M) ^ U 
dn 2 dn 2 



(4.14) 



For uj, M <C 1 the solution of this continuum equa- 
tion is e xpe cted to coincide with the exact eigenfunction 
of Eq. (4.6) for n 3> 1. In Appendix B, we verify this 
for the special case M = 0, where it is possible to solve 
directly the discrete Schrodinger equation. 



The continuum differential equation (4.14) must be 
supplemented by an appropriate boundary condition. A 
natural physical requirement is that <p{n, M) — > for 
n — > oo. Because we have still have the freedom to choose 



ci and C2 in Eq. (4.13), the normalization at the origin 
can remain arbitrary at this stage. 

The two constants in Eq. ( 4.13| ) then require tw o ad - 
ditional constraints. The first comes from Eq. (4.9), 
<h = 1, giving 



(4.15) 



ci^(l,M) + C30(l,-M) = 1. 



The second constraint is obtained from Eq. ( }4.6| ) for 
n = 1: 



-M(<fi 2 - <t>o) +U(f>l 



(4.16) 



In the limit M, w<l, and using Eq. ( 4.13 ), this becomes 
c 1 ^(l,M) = c 2 0'(l,-M). (4.17) 



D. Solution 

We are now in a position to obtain the solution to the 
continuum equation. Under an exponential cha nge o f 
variables, n = e z , with $(z,M) = <j)(e z ,M), Eq. ( 4.14| ) 
takes a more illuminating form, 



dz 2 



M- 



dz 



0. 



(4.18) 



Indeed, for M = this is equivalent to a Schrodinger 
equation for a particle moving in an exponential po- 
tential. Since z = Inn, the domain of the equation is 
< z < oo. For small u>, however, the potential is neg- 
ligible for small z. It rises very abruptly and becomes of 
order one for z ~ z,„, where 



= line 



(4.19) 



In the small uj limit of interest, th en, th ere is a region 
of divergent width over which Eq. ( 4.18 ) effectively de- 
scribes a free particle. The width of the potential is, 
however, only of order one. Thus we expect that on the 
scale of z w , the exponential potential can be well approx- 
imated by a hard wall. 

In the following we make this hard- wall approximation, 
replacing the continuum equation ( [4.18 ) by 



— -M— 

dz 2 dz 



$ = 0, 



(4.20) 



with < z < z w . and the boundary condition 

*(*«,) = 0. (4.21) 

This simple equation can be readily solved (below) and 
the density of states extracted. 

As shown in Appendix C, the hard-wall approximation 
is not necessary, since the exact continuum equations can 
be solved explicitly. Although the hard-wall wave func- 
tion differs from the exact solution, the resulting density 
of states coincides, up to an overall (non-universal) mul- 
tiplicative constant. 



The solution of Eq. ( 4.20) con sistent with the hard- wall 
boundary condition (Eq. (|4.2l|)) is 



Mz v 



V 



(4.22) 



where we have assumed that ui,M <C 1, and chosen a 
convenient (but arbitrary) normalization. As M — > 
this reduces to 



$0) =1 . 



(4.23) 



Applying the constraints in Eqs. (4.15,4.17), determines 
the constants as 



ci = c 2 



z,„Muj m 



1-LU 



2M 



(4.24) 



To extract the density of states, we use Eq. (4.11) to 
write, 

/>oo 

G{0, ioj) = 2iac 2 / dn (f>(n, M)<j>(n, -M) 



2i Cl c 2 I dze z $(z,M)$(z,-M). (4.25) 
lo 
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The integral can be readily performed giving, 
G(0,iu)=i—f(u, M ), 



with the exact scaling function 

Y 



f(Y) 



1-Y 2 



(4.26) 



(4.27) 



Our present implementation of the hard-wall approxima- 
tion reproduces the exact value of the non-universal con- 
stant, c = 4, obtained in Appendix C, but this result 
depends upon the precise position of the wall. 

An exact expression for the density of sta tes in the crit- 
ical regime can now be obtained from Eq. ( |2 . 25| ) by per- 
forming an analytic continuation, p(e) — (l/Tr)ImQ(0; e+ 
i0 + ). Noting that Q(0,iuj) is pure imaginary and odd in 
to, one readily obtains (for e > 0), 



with 



P (e) = —g(e^), 



1 Y 2 (l + Y 2 ) 

fP0 = ^ = ^ 



In the M — > limit, this reduces to, 



I In el 3 ' 



(4.28) 



(4.29) 



(4.30) 



a result obtained previously by other methods. 

A special feature of Id, is that the typical localiza- 
tion length, £, can be extracted from the real part of the 
Green's function at x — 0. As derived in Ref. [24|, £(e) 
satisfies 



dt 1 



P / de 



-ReQ{x = 0;e) 

7T 



(4.31) 



Performi ng a n analytic continuation to real energy us- 
ing Eq. ( fl.26|) gives ReQ(0;e) ~ — 1/cj lne| 2 . Integration 
from e = using the fact that £ -1 (0) = gives the result 
£~|lne|. 



V. FERMION GREEN'S FUNCTION 



To determine the mean correlation length, we need to 
calculate t he G reen's function, G(x, iui) at non-zero x. 
From Eq. (3.23) this takes the form 



G af} {x,iuj) = ^(-l) Q L (0| F a e- xIi Fl\Q) R . (5.1) 

Given the quantum numbers of the ground state, the 
state F^\0) H clearly has N B = -1, N F = 2, T = N = 1, 
r = 0. We can therefore insert the resolution of the 



identity from Eq. (3.81) to give 



G a0 (x,iLo) = i(-l) a Y, l(0\ F a \ - 1,2,1,0, E) R L { ~ l,2,l,0,£|fj|0) fl e- 



(5.2) 



E 



Our task is thus to determine the matrix elements and eigenvalues needed to carry out this sum. To do so, we expand 
the eigenstate in the appropriate basis, 



- 1,2,1,0, E) R = J2xn\n) 1 , 



n=0 

where we have abbreviated 

|n)i = |- 1,2, 1,0, 1/2 + n). 
The wavefunction Xn then satisfies the Schrodinger wave equation, 



2u> 



Xn+2M [nxn-i - (n + l)Xn+l] 

-(n + 2)(n + l)Xn+2 + (2n 2 + 2n + l) Xn - n{n - l) X n-2 = E X , 



(5.3) 



(5.4) 



(5.5) 



where we have again set g = 1. Again, this equation has 
the property that multiplication of a solution by (— 1)" 
yields a solution for M — > —M. We therefore expect 

Xn(M,E)=c 3 x(n,M,E) + Ci (-l) n x(n,-M,E). (5.6) 



2n- 



dn 



2M 2n 



dn 



Naively, x( n i M, E) can be obtained in the continuum 
limit. Converting finite differences to derivatives, we ob- 
tain 



2un 



X (n,M,E) = E X (n,M,E). 



(5.7) 
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Upon tranforming to logarithmic variables, 

x(n,M) = $ E (z = hm,M), 

one has 

(2d z + i f - 2M (2d z + 1) + 2uje z ] <S> E (z, M) = E<S> E (z, M) 



(5.8) 



(5.9) 



Based on previous experience, we expect that the 
hard-wall approximation gives exact results in the scal- 
ing limit. In Appendix D we verify this explicitly, by 
cons tructing exact solutions of the continuum equation 
(5.7). Within the h ard- wall approximation the potential 
2uie z — > in Eq. (5.9) is dropped, and replaced by a 
boundary condition at z w = | kiwi, 



$ E (z w ,M) = 0. 



(5.10) 



The general solution of Eq. (5J)) with the hard- wall 
boundary condition is 



* B (z) = e"( 1+M >/ 2 sin (P(z - z w )/2) . 



(5-11) 



where (3 = y/E - M 2 . 

So far, we have not determined the spectrum , or al- 
lowed values of f3. In an ordinary quantum problem, 
these eigenvalues would be fixed by a boundary condition 
at z — 0. In this case, however, such a simple treatment 
is problematic. The difficulty arises because, unlike in 
the ground state sector, neither \ n = 1 nor Xn = (— 1)" 
are solutions in the limit ui = M = 0. We therefore ex- 
pect a non-trivial solution for n = 0(1), in which the 



discreteness of n is important and the continuum limit is 
not valid. 

Fortunately, for 1 ^ n « we can obtain an 

asymptotic approximation which does not rely upon the 
continuum limit. This is possible beca use fo r n -C l/u>, 
the loJ z term in the Hamiltonian Eq. ( |3.31| ) can be re- 
garded as a small perturbation. Neglecting this term, H 
is a function only of J x , so that the eigenfunctions of H 
are simply eigenfunctions of J x . As a first step, consider 
the state \a) with 



J x \a) = ia\a). 
Expanding \a) in the J z basis, 



= ^2^n{a)\n) l , 



(5.12) 



(5.13) 



one finds the simpler Schrodinger equation 

(n + 1)^71+1 (a) - mj} n -x{a) = -2iaip n (a). (5.14) 

This equation can be solved exactly (see Appendix E). 
For large n, the solution which is well-behaved at the 
origin behaves asyptotically as 



-1/2 



[(2n)- M r(l/2 + to) + (-l)"(2n) m r(l/2 - ia)] 



(5.15) 



Eigenstates of H therefore take the form 

Xn ~ c+ip n (a + ) + c-ip n (a-), for 1 < n < 1/oj, 

(5.16) 

where a± are the two roots of the equation 4a 2 +4iMa = 
E, i.e. 



a± 



-iM±/3 



Xn - e- z ' 2 [e- Mz ' 2 



(5.17) 

c+e -0z/2 + c _^z/2 + (_ 1) n e M,/2 



and 13 = y/E- M 2 . 

Mathematically, Eq. ( 5.16| ) is an outer solution, valid 
outside a boundary layer which occurs for large n. To 
obtain a complete solution, it must be matched to the 
inner solution, valid "inside" the boundary layer, which 
is just the continuum regime of large n. Within the 
hard-w all ap proximation, this is just the standing-wave 
in Eq. ( pTl ). 

To match the two solutions, we let n — e z in 
Eqs. ( |5~T5| , p^ ), which gives 



c+e 



ipz/2 



+ c_e 



-if3z/2 



Similarly, using Eqs. (5.6,5.11), the continuum solution gives 

Xn - e~ z ' 2 {c 3 e- M2 / 2 + (-l)"c 4 e Mz/2 } sin (fi(z - z w )/2) . 



(5.18) 



(5.19) 
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These expressions are equal in two situations. One can C3 = (— 1) + C4 = c. 

take c+ = c_ =03/2 = 04/2 if /Jz^, is an odd multiple of 

7T. Alternatively, c + = -c_ = c 3 /2 = -c 4 /2 if /3z w is an The discrete quantum number fc = 1, 2, ...00 
even multiple of tt. The final, matched solution for both 
cases can thus be written 



(5.24) 



*<*>(*, M) = <S> Ek (z, M) = e -( 1+M >' 2 sin {faz/2 + k ) , 

(5.20) 



The correspond ing l eft eigenstate is obtain by acting 
with U from Eq. fl3.74j ). Defining 



where 



\k) R/L = \-l,2,l,0,E k ) R/L , 



we choose 



fi k = Ttk/z w , 



E k = M + — 



7fc = 



k even 
tt/2 k odd 



(5.21) 
(5.22) 

(5.23) 



\k) L = (-l)^e 



fc+l iir(J z -l/2) 



\k) R . 



(5.25) 



(5.26) 



With this choice, normalization implies that the constant 

— 1/2 

c = z w ■ Thus the final expressions for the (normal- 
ized) eigenstates are 



\k)R = ^ 1/2 E [*(™' M ) + (-l)" +fe+1 x(",-Af)] In)!, 

n 

\k) L = z- 1 / 2 Y / [(-V n+k+1 x(n,M)+ X (n 1 -M)] |n) x . 



(5.27) 
(5.28) 



These actually form an orthonormal set 
L ( k\k') R = S k k', 



as can be verified by direct computation. 



(5.29) 



Having obtained the full set of eigenvalues and left and righ t eigenfunctions of H in the appropriate sector, it is a 
simple matter to evaluate the Green's function using Eq. (5.2). Using the hard- wall eigenfunctions, we find 



G a p{x, iu>) = i 



32tt 2 



M 2 lu 2M 



i-ir 



uj\\noj\ 3 (1~lu 2M ) 2 

OO 

< ^(-l) fc+1 fc 2 \(-l) a u; M / 2 + (-l) fc+1 ^- M / 2 l \{-lfw M ' 2 + (-l) k+1 L 



-M/2 



-Tr 2 k 2 x/\ lnw| 2 



(5.30) 



k=l 



One thereby obtains the final form for the Green's 
function, exact in the scaling limit: 



Q(x, iu) 



even 



(5.31) 



where the universal scaling functions are given by, 



£00 = YfZ(X) = 



YlnY 



l-Y 2 
F*(Y) = Aj2k 2 o- k2Y 



k=l 



F°{Y) = A^{-l) k+1 k 2 e- k 



-krY 



(5.32) 
(5.33) 

(5.34) 



k=l 



Here the non-universal amplitude A = 32tt 2 within the 
hard-wall approximation, and we have defined two cor- 
relation lengths, 



l/M 2 



L. = [ — 

7T 



(5.35) 



Notice that Q(x;iu>) is pure imaginary for x even and 
pure real for x odd, as dictated by particle-hole symme- 
try. Eq. (5.31) is valid for lu > 0. Particle-hole symmetry 
then determines the Green's function for u> < 0, since 
G(x\ itu) is odd in u> for even x and even in u for odd x. 



Eq. ( [5.35j ) gives us the correlation length exponent, 
v = 2, defined by £ ~ M~", and the dyna mical expo- 
nent, z = 00 , defined by lu ~ £ z . Eq. (5.31) is actually 
simpler than one would generally expect on the basis of 
scaling - the most general scaling form would not factor- 
ize as it does here. 

The correlator C(e) = (l/n)lmQ(e + i0 + ) can be ob- 
tained via analytic continuation, using the symmetry 
properties of the Green's function under lu — > — to. One 
finds 
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~ x/iM , (5.36) 



(_l)n+l f x 

C(x = 2n + 1, e, M) = K . ,. F° -, e 



p -x/Sm ) (5.37) 



where the energy dependent "localization" length is 



6 = 



In i 



(5.38) 



The even sublattice scaling function no longer factors, 



F*(Y,Z) = f°(Z) 
1 



-F-(Y)+YF-'(Y) 



Z]nZf*'(Z)FZ(Y), 



(5.39) 




while the odd sublattice scaling function remains simple, 



F e °(Y,Z) = FZ(Y)K(Z). 



(5.40) 



Fig. 2: Scaling functions for the Fermion Green's func- 
tion. The solid and dashed lines are the scaling functions, 
F e e (Y") and F°(Y), for correlators C(x) between two sites 
on the same and different sublattices, respectively. 

The different sublattice Green's function behaves quite 
differently at short distances. Using the Poisson summa- 
tion formula, it is straightforward to show that 



The scaling forms, Eqs. ( 5.36 - 5^401 ), encode several sig- 
nificant physical properties. First consider the same sub- 
lattice correlation (x = 2n), for simplicity at zero stag- 
gering (M — 0). For distances shorter than the cor- 
relation length, this has a slow power-law decay, since 
F e e (F, 1) - Y- 3 / 2 , for Y < 1. In particular, 



C(x = 2n, e, M = 0) 



(-1)* 



e|lne| 3 \x\ 3 / 2 '' 



(5.41) 



2 y3/2 



7T 

1 

4Y 



(5.42) 



for Y 1. This implies that C(x = 2n + 1) is much 
smaller than the same sublattice correlator . The sup- 
pression is due to the fact that sublattice mixing of wave- 
functions occurs only at finite e. The exact zero energy 
"wavef unction" , discussed in detail in the next section, 
lies entirely on a single sublattice. 



for | a; | <C £ c . This can be understood as the product 
of the density of states and a two-point "wavefunction 
correlation" , with multifractal scaling exponent (see the 
next section) y(q = 1) = 3/2. For distances x » £ c 
(and M = 0), even the rare wavefunctions are localized, 
and C{x) decays exponentially. The full scaling function, 
which describes the crossover between these two limits, 
is plotted in Fig. 2. Note the change of sign for Y ~ 2.5 
- this may be interpreted physically as arising from the 
first node in the dominant rare wavefunction at energy e. 



VI. MULTIFRACTAL WAVEFUNCTIONS 

The above results reveal that the ensemble averaged 
single particle Green's function decays as a power law 
with an exponent 3/2 for distances smaller than the 
correlation lenth. When m = 0, the relevant scale is 
£ e = | lne| 2 /7r 2 , which diverges at the band center. Such 
power law spatial scaling is in striking contrast to the 
behavior of the typical Green's function, which decays as 
a stretched exponential, Gt yp . ~ exp(— ax 1 / 2 ) at critical- 
ity. This stretched exponential form follows directly from 
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the exact zero energy wave function, Eq. (2.13), whose 
logarithm undergoes a one-dimensional random walk. 

Power law scaling of the average Green's function at 
criticality can al so be understood (albeit more subtly) 
in terms of Eq. (2.13). To see this it is useful to con- 
sider ensemble averaged correlation functions of the zero 
energy wavefunction. As with higher dimensional local- 
ization transitions, such as the plateau transition in the 
IQHE, this critical wavefunction is expected to exhibit 
multi-fractal scaling characteristics. As shown in a recent 
paper, such correlators in this Id case can be computed 
exactly, via a mapping to Liouville quantum mechanics. 
This mapping exploits the equivalence between imagi- 
nary time quantum mechanics and the one-dimensional 
random walk. Slightly generalizing this work, we com- 
pute below the full multi-fractal spectrum for the Id crit- 
ical case. This calculation is instructive since it reveals a 
link between the supersymmetry calculations and Liou- 
ville quantum mechanics. 

To extract average wavefunction correlators, it is nec- 
essary to consider normalized states. We thus consider 
a finite system of len gth L, and normalize the wavefunc- 
tion from Eq. ( |2.13| ) over the finite interval |x| < L/2. 
Focussing on one component of the spinor wavefunction, 
say <j) + (x), an appropriately normalized wavefunction, 
denoted as ip(x), can be written, 



with normalization 



r L/2 

Af = / dxe 2z ^ x \ 

-L/2 



(6.1) 



(6.2) 



where we have defined d x z{x) = m(x), with the random 
potential m(x) assumed as before to be Gaussian with 
[m(x)m(x')] ens , = 2g8{x — x r ). It remains to specify the 
boundary conditions on ip(x) at x — ±L/2. For technical 
reasons it is convenient to "pin" the Id random walker 
at the ends, taking z(x = ±L/2) = 0. 

We focus on the ensemble averaged correlation func- 
tion, 



W q (x ) L) = Mx)m\ q ] 



(6.3) 



between two points separated by a distance x, assumed 
much smaller than L. The one-point function W q (L) = 
yV q {x = 0, L), referred to as a participation ratio, is ex- 
pected to vary as a power law with system size: 



W q (L) 



1 



(6.4) 



with d = 1 the spatial dimensionality. The exponent 
r(q) is often written as r(q) = (q — l)D(q). For a plane 
wave or non-fractal wavefunction (as in a 3d metal, say) 
D(q) = d. For an exponentially localized wavefunction, 
D(q) = r(q) = 0. A simple fractal would be charac- 
terized by a g-independent D, different that the spatial 



dimension, whereas in a multi-fractal D depends on q 
and equivalently r(q) is a non-linear function of q. 

The two-point function W q {x 1 L) for x much larger 
than microscopic lengths (ie. lattice spacing) yet much 
smaller than L is also expected to exhibit power law scal- 



ing: 



W g (x,L) 



1 



ljd+T( q ) x y( q ) ' 



(6.5) 



with a spectrum of exponents y(q). For some multi- 
fractals a relation can be obtained between r(q) and y(q), 
but in general they can be independent exponents. 

To extract r(q) and y(q) for the Id critical wavefunc- 
tion, we follow closely Shelton and Tsvelik.E3 The corre- 
lation function can be expressed as a functional integral 
over the random walk configurations z(x) as, 



with 



W q (x,L) = ±- I Vz\ij(x)i;(0)\ie- s °, (6.6) 



Z = I Vz e~ s ° 



and an action, 



i r L/2 

S = — dx {d x zf 

4 5 J -L/2 



(6.7) 



(6.8) 



Here the prime on the integration indicates the bound- 
ary conditions, z(x = ±i/2) = 0. In the following we 
put g = 1. The functional integral over disorder con- 
figurations is non-trivial due to the normalization of the 
wavefunctions. Fortunately, the normalization can be ex- 
ponentiated via the identity, 



1 f°° 

r w Jo 



(6.9) 



where it can be absorbed into the action. In this way one 
obtains, 

1 C°° 

W q (x,L) = — <W-V* (a V z(0) ), (6.10) 

r(«) Jo 



where the average is given by 

^qz(x) e qz(0)^ = _L ( Vz & qz(x) e qz(0) 



Zo 



(6.11) 



with the total action 

rL/2 
-L/2 



r L/2 i 

S= / dx[-(d x zf +toe 2z ]. (6.12) 

J-L/2 4 



If x is viewed as an imaginary time coordinate, this 
average is seen to be equivalent to a path integral rep- 
resentation of the quantum mechanics of a particle with 
coordinate z moving in an exponential potential. Passing 
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to the operator representation of this quantum mechanics 
by defining a quantum Hamiltonian, 

H = -d 2 +u;e 2z , (6.13) 

the above average can be written as a quantum expecta- 
tion value, 



-LH/2 e qz(x) e qz e -LH/2 



|0) 



<0|e-"*°|0) 



.14) 



with e qz ^ 



e xH e qz t 



-xH 



and Hq = —8%. Here |0) is a 



position ket with z = 0. 

Evaluating the wavefunction correlator has thus been 
reduced to solving for the quantum mechanics of a parti- 
cle moving in an exponential potential - Liouvillc quan- 
tum mechanics. This form is identical to that which arose 
in the bosonic sector of the supersymmetric calculation 
of the Fermion Green's function. In that case, the co- 
ordinate z was related to the boson number via n = e z , 
and u) was a small imaginary part of the energy. The 
supersymmetric calculation can thus be viewed as a su- 
persymmetric version of Liouville quantum mechanics. 
In earlier work, Kogan, Mudry and Tsvelik showed that 
the wave function correlators for a two-dimensional par- 
ticle decribed by a Dirac equation with random vector 
potential (for which the exact zero energy wavefunction 
can also be written down explicitly) could be formulated 
in terms of Liouville field theory. Perhaps such a 2d lo- 
calization critical point can be formulated in terms of a 
supersymmetric Liouville field theory. 

For simplicity we evaluate the above quantum expec- 
tation value within the hard-wall approximation, which 
should give the correct scaling behavior for the wavefunc- 
tion correlator. As before, we replace the exponential 
potential by a hard- wall at z w , with we 2z " = c, for a 
constant c of order one. The value of c affects the overall 
prefactor in the correlator. We choose c = 1/2 which 
gives the correct normalization, W q= i(L) = 1/L. Since 
the quantum particle is contrained to have z < z w , when 
uj > 1 (and z w is negative) the particle's wavefunction 
vanishes at z = 0, so that the quantum expectation value 
in Eq. (|6.14 ) vanishes identically. We can thus restrict 
the integration over uj to uj < 1. 

The denominator in Eq. (|6.14 ) is the propagator for a 
free random walker (i.e without the hard-wall) and can 
be readily evaluated giving (47rl/)" 1/2 . To evaluate the 
numerator it is convenient to let Z <* ZfjQ z ^ SO that 
the quantum particle is then constrained to have z > 0. 
Within the hard-wall approximation, the correlator can 
then be expressed in terms of the free Hamiltonian as 



W q (x,L) 



V / 4ttX 



, 1 ^<*J e -^o/2 e -^x) e - g * e -W2u \ 



(6.15) 



with z restricted positive and z w = |lnw|/2 (the factor 
of two difference between this definition and the one used 



in the SUSY calculations is a consequence of a different 
choice of normalization of the field z(x)). To evaluate this 
quantum average we introduce a complete set of standing 
waves 



(6.16) 



(z\k) = \-sm(kz), (k>0) 



which are eigenstates Ho\k) = k 2 \k), and appropriately 
normalized on the interval z > 0: (k\k') — S(k — k'). 
Inserting the resolution of the identity, 



1 = / dk\k){k\, 
'o 



(6.17) 



into Eq. (6.15), evaluating the matrix elements in closed 
form and performing the integration over u> gives for large 
L the final result: 



W(q 2 x) 1 

29-vr(g) V 



with 



W(x) 
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dkk 2 



(1 



(6.18) 



(6.19) 



The crossover function W(x) interpolates between one 
at x — and W(x) ~ x~ 3 ^ 2 for x » 1. The crossover 
length, x c = 1/q 2 , is the characteristic distance over 
which the logarithm of the wave function ip(x) changes by 
order 1/q (or equivalently the "time" it takes the Id ran- 
dom walk to move a distance Sz ~ 1/q)- This crossover 
scale clearly depends on the strength of the disorder. For 
the original lattice tight binding model, when the random 
hopping strengths are comparable to the mean hopping 
strength, St n /t of order one, this crossover length is of 
the order of the tight binding lattice spacing. Thus it is 
clear that the form of the crossover function W(x) for 
x of order one, cannot be universal. In fact, the precise 
form above is particular to the hard wall approximation, 
and an evaluation using the exact Liouvillc eignenfunc- 
tions would give another form, although they agree in 
their universal large x behavior. 

For large x, W g (cc, L) ~ x~ 3 ^L~ x . By comparing with 
Eq. ( |6.5D , the exponents r(q) and y(q) are seen to be in- 
dependent of g, with r(q) = and y(q) = 3/2. A van- 
ishing r(q) is characteristic of an exponentially localized 
wavefunction. However, the two-point correlator for an 
exponentially localized wavefunction also decays expo- 
nentially (y(q) = oo), in contrast to the present Id wave- 
function which exhibits power law correlations. The Id 
critical wavefunction is typically quasi-localized (centered 
around a maximum) with stretched exponential decay. 
However, the average two-point correlator at separation 
x is dominated by the rare wavefunction which has a sec- 
ondary maximum close in magnitude but separated spa- 
tially (by distance x) from its primary maximum. The 
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likelihood of this involves the extremal statistics of a Id 
random walk near a global maximum (absorbing wall), 
which is being described mathematically above by quan- 
tum mechanics near an exponential (or hard wall) poten- 
tial. As expected, the two-point wavefunction correlator 
at q = 1 decays with the same exponent, y(q = 1) = 3/2, 
as the average Green's function obtained with supersym- 
metry. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we have presented a detailed SUSY anal- 
ysis of the critical properties of the zero energy dereal- 
ization transition in the Id random hopping model. This 
Id random critical point has been of interest for many 
years, originating with the pioneering paper by Dyson 
in 1953 on a related model of a Id harmonic chain with 
random spring constants.E3 Most of the prior work, using 
a variety of different approaches, has focussed on prop- 
erties derivable from the mean local Green's function, 
specifically the density of states and the typical local- 
ization length, By employing a novel real-space RG 
approach to analyze phase transtions in a class of closely 
related random spin-chains, D.S. Fisher has recently ex- 
tended the analysis to extraet-the spatial dependence of 
mean correlation functions.r 9 lr°l An important element is 
the emergence of a second correlation length, £, which de- 
termines the spatial decay of mean (rather than typical) 
correlation functions. To the growing body of knowledge 
concerning this Id random critical point, we add several 
new results in this paper, (i) Using SUSY we have com- 
puted the exact two-parameter scaling functions for the 
mean Fermion Green's function, (ii) By employing Li- 
ouville quantum mechanics, we have extracted the set of 
multifractal scaling exponents r(g) and y(q) which char- 
acterize the critical wavefunction pair correlators. 

Together, these two results encapsulate the important 
universal scaling characteristics of this Id random criti- 
cal point. The spatial dependence of the mean Fermion 
Green's function is controlled by two lengths, a mean lo- 
calization length which diverges upon approaching the 
band center as £ £ ~ |lne| 2 , and a mean "staggering 
length", varying as £m ~ M~ 2 when the strength, M, 
of a staggering in the hopping strengths is taken to zero. 
These two lengths are to be contrasted with their coun- 
terparts, denoted £ e and £m, which charaterize the spa- 
tial decay of the typical (rather than ensemble-averaged) 
Greens function. From the singular behavior of the den- 
sity of states, one can infer that the typical localization 
length diverges more slowly, as £ e ~ \lne\. Likewise, the 
typical staggering length which follows rather directly 
from the nature of the exact (decaying) zero energy wave- 
function, diverges more slowly, £m ~ M _1 , than it's 
mean counterpart. 

For spatial separations x <C £ e , £m, between two points 
on the same sublattice, the mean Fermion Green's func- 



tion varies as an inverse power law of x with universal 
exponent 3/2. This result also follows from an analysis of 
the zero energy wavefunction whose logarithm undergoes 
a Id random walk, with the exponent 3/2 being related 
to extremal properties of the random walker. In contrast, 
the typical Green's function for x <C £ £ ,£m is expected to 
decay as a stretched exponential, Gtyp ~ e~ c ^ , reflecting 
the typical behavior of the random walker. 

A key motivation for the present paper was to in- 
vestigate in detail the novel features which emerge in 
a SUSY formulation of a random critical point. The 
calculation proceeded by expressing mean correlators in 
terms of quantum mechanical expectation values for a 
SUSY Hamiltonian, which involved a single superspin. 
This Hamiltonian has a number of notable features: (i) 
It is non-Hermitian, with distinct left and right eigen- 
functions. (ii) It has a unique zero energy ground state, 
as dictated by supersymmetry, and the excited states are 
organized into supersymmetric doublets and quadruplets, 
(iii) The right (or left) eigenstates alone do not span the 
Hilbert space - the Hamiltonian is thus "defective" . (iv) 
The Hilbert space is infinite, due to the non-compact 
SU(1,1) bosonic subalgebra of the superspin group, (v) 
The eigenstates explore the outer reaches of the non- 
compact manifold, in a manner which can be described 
by Liouville quantum mechanics. 

It is our hope that a thorough undertanding of 
these unusual features will be helpful in extending the 
SUSY approach to attack two-dimensional random crit- 
ical points, such as the IQHE plateau transition. It is 
tantalizing to speculate that some appropriate supersym- 
metric version of Liouville field theory might give a cor- 
rect description of derealization transitions in 2d. 
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APPENDIX A 

The mapping from the random transverse field Ising 
chain to a free Fermion model was introduced by Shankar 
and Murthy. We briefly recapitulate this mapping. In 
terms of Majorana Fermions, 

= :ti n « t 7 - 1 ) 

* m<n 

* m<n 

which satisfy fa, m , ru n } — Sij5 mn , the random Ising 
Hamiltonian, Eq. J2.27| ) can be re-written as, 
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2iK l!n r]i >n ri2,ri + 2iK 2 .n'lll,n' l l2,n+l 



(7.3) 



A continuum limit can be taken by putting x = ndx, 
K\ = dx/2 and K2 — (1/2 + m(x))dx, and converting 
the sums to integrals. This gives 



dxrj 



a x id x + m(x)a y 



V, 



(7.4) 



where we have defined a two-component Majorana field, 
V = (771 7 772)- For spatially uniform m this model de- 
scribes criticality in the pure 2d Ising model, with the 
phase transition occuring at m = 0. 

To complete the mapping it is convenient to consider 
a path integral representation of the partition function, 
Z = Trexp(—f3H), which can be written as a functional 
integral over Grassmann fields, t)(x,t), with associated 
Euclidian action: 



S = J dx^-n(x,uj) 



a x id x 



ma" 



r)(x,-w). (7.5) 



These can be decomposed into new Grassmann fields by 
defining, 



Va (w) = i> a (OJ) , r) a (-w) = 1p a (w), 



(7.6) 



for positive lu and a = 1, 2. In terms of these new fields 
the action becomes, 



S- I ^S u , 



with 



dxip 



o~ x id x + m(x)a v + iuj 



(7.7) 



(7. 



Notice that the functional integral factorizes into a prod- 
uct over independent frequencies. In the following we 
focus on only a single frequency. The action at a single 
frequency can be cast into the form of Eq. (2.6) by a ro- 
tation in "spin-space" around the y-axis by 7r/2, which 
takes o~ x — > a z , giving 



S' 



dxip(h + iiS)tp, 



(7.9) 



with h the Id random Hamiltonian in Eq. (2.7). Note 
that, in this case, a non-zero mass mo corresponds sim- 
ply to the deviation from the Ising critical point. 



APPENDIX B 

For the special case of the ground state wavefunction, 
we are in fact able to obtain an exact solution without 
taking the continuum lim it. T his is possible because the 
difference equation, Eq. (|4.6|), is linear in n. Here, we 



specialize to the case n = 0, in which this solution is 
especially simple. 

For M = 0, the Schrodinger equation decouples on 
even and odd sublattices and can be solved independently 
on each. To bring this out, we define 



for n = 0,1,2,. 
then obey 



7n: <hn+l = 7«, (7-10) 

.. The even and odd sublattice fields 



0,7° = (n + 1) 7 ° +1 - (2n + 1) 7 ° + n 7 °_i n > 0, 
uy e n = (n + l/2) 7 ^ +1 - + (n - l/2)-y*_ 1 n > 0. 

To solve them, we define the generating function 

W) N . (-El ) 



where P — o, e. Consider first the odd sector. Multiply- 
ing the equation for 7° by w n and summing gives 



(1 - w + w) 7 ° 



(7.12) 



This is easily solved by separation of variables, to give 
C 



7 



1 



■ exp 



1 



(7.13) 



where C is an arbitrary constant. Not the strong diver- 
gence at w = 1. This implies unacceptable behavior for 
7° at large n. 

The even sector is (fortunately!) rather more co mpli - 
cated. The crucial difference is the fact that Eq. (4.6) 
is valid only for n > 0, leaving an extra free parameter. 
Carrying out the transform in this case gives 



2 ^7 e 











R) 


— UJ 









1 / 1 



2U" 71 



(7.14) 



Where we have imposed the normalization 7 § = 1. Note 
the appearance of 7' as a parameter in the equation. It 
must be adjusted to achieve a well-behaved solution. 

This inhomogenous equation can be solved by intro- 
ducing the integrating factor (1— w)/y/w exp[— u>/ (1— id)]. 
The solution is 



1 



dy 

y/vQ- - y) 



2y 



e -«/(i-») ) ( 715 ) 



where 8 = 1 — 7 |. Performing an integration by parts 
leads to the form 



1 



w 



1 + 



J{W) 



(7.16) 
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where 
J{w) = 



dy 



o y/vO- - v) 



1-y 2 

(7.17) 

To avoid the strong divergence of 7^ , we clearly need 



J(1) = 0. 

This fixes 5. This implies that 



(7.18) 



J(w) = - 



dy 



Vtf(l-V) 



l-y 2 



S 

LU 



-U/(1-V) 



(7.19) 



We now change variables via y = 1 — u>/t, and also define 
s = (1 — w)/lu. Then 



J = — 



dt 



/a ty/1 — w/t 



t U) 

2 



(7.20) 



We may now take the limit lu — > 0, with s fixed, and 
S > uj. Then 



where 



E x {x) = I dte-t/t 



(7.21) 



(7.22) 



is the exponential-integral function. Plugging back in 
gives, finally, 

rW _ 5 _«V. ft( l / .,__jf (7.23) 

In fact, 7 e (s) is nothing but the Laplace transform in the 
limit lu -C 1, 



7 [s 



E 



e —nus 

7„e 



dn7 e (n)e-"" s = [L 7 e ](ws). (7.24) 
We can therefore invert it using the inversion formula 

f c+zoo 



7n 



y . 
- y . 



5 
2 Jo 



die"' 



1 



J uis i 



27TI 1 + ts 



2./0 "t 

2 Jo 



dt-e-™ 1 *- 1 



-nt—uj jt 



(7.25) 



This is precisely the modified Bessel function solution 
obtained from the continuum limit (see appendix E). 



APPENDIX C 

In this Appendix we obtain the density of states with- 
out resorting to the hard- wall approxima tion, by solving 
exactly the full continuum equation Eq. (4.14). This can 
be accomplished by employing an inverse Laplace trans- 
form, defining 



(j)(n,M) 



dte~ nt j>(t). 



(7.26) 



'0 

Provided t 2 (p(t) — > as t — ► (as required for a well- 
behaved solution as n — > 00), insertion into Eq. (4.14) 
leads to the simple transformed form 



.2 # 



1 ^ + < 1 



This has the general solution 

m " 



j in 



l-M 



(7.27) 



(7.28) 



where a is an arbitrary constant. The "unnormalized" 
wavefunction is thus given by 



4>(n, M) = a 



dt 



^ — nt—uj jit 



(7.29) 



which is the integral representation of a Bessel func- 
tion, 4>{n,M) = 2a(uj/2n) M / 2 K M (\ / 2u~m)- To evaluate 
Eqs. (4.15|,|4.17[), we need 0(1) and </>'(!). These are de- 



termined by making the change of variables x 
which yields 



-M 



dx exp 



(xW~ 



-l/M 



(7.30) 

with the scaling variable W = (tu/2) M . In the scaling 
limit lu,M — > with W fixed and finite, Eq. (7.30) can 



be simply evaluated. Since each of the arguments in the 
exponential goes to zero or infinity, the limits of integra- 
tion are restricted giving, 



4>{n 



M 



(7.31) 



(using 2 M w 1). The same change of variables can be 
used to extract the n derivative, giving in the scaling 
limit 



dn 



n=i = -a. 



(7.32) 



Comparison with the hard-wall forms shows that the 
constants c\ and C2 are identical provided we take a — 
l/z w = \ lncj| _1 . 

With the exact solutions of the continuum equations 
in hand, one can readily evaluate the densit y of states 
by inserting the integral representations Eq. ( 7. 29] ), into 
the expression for G{ito) in Eq. (4.25). The n-integration 
can be readily performed. In the scaling limit M <C 1 
the remaining two ^-inte grations a re simple and yield an 
identical result to Eqs. (4j§-fL27|). 
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APPENDIX D 

In this appendix, we obtain the exact excited state 
wavcfunctions in the continuum limt in the Np = 0, 2, 
Nb = — 1 sectors, and show that they lead to the same 
scaling form for the Green's function as does the hard- 
wall approximation. Beginning with Eq. (5.7), we make 
the change of variables 



x(n,M) = a(l3)n-( 1+M >/ 2 x(n,(3), 



(7.33) 



where a(/3) is a normalization constant to be chosen later 
in order to maintain the closest possible agreement with 
the hard-wall solutions in section V. The transformed 
wavefunction then satisfies the simpler equation 



1 



dn 2 dn 



+^ + i(/ ?2 - 2 H 



(7.34) 



where f3 = VE - M 2 . Eq. (Ol is a standard equa- 
tion of classical mathematical physics. Its solutions are 
modified Bessel functions of imaginary index: 



(7.35) 



where we have chosen the solution K which decays at 
infinity. Note that we have assumed E > M 2 , for which 
(3 is real. It is straightforward to show that there are no 
satisfactory solutions with E < M 2 . Very few results are 
readily available for these functions at imaginary index. 
Instead, we will make heavy use of the integral represen- 
tation, 



la 



/2nu cosh t 



cos f3t, 



(7.36) 



whic h ca n be verified by direct substitution into 
Eq. ( 7.34 ). A se cond useful form is obtained by inte- 
grating Eq. (7.36) by parts: 



\Z2nw 



— y2nuj cosh t 



sinhisin/3i. (7.37) 



P Jo 

Performing the integral over n gives 

2 f°° ,., . „. . „., sinhisinhi' 



In 



1 



dtdt' sin/3tsin/3t' ■ , , 

o (cosht + cosht')^ 



The first task at hand is to determine the spectrum, or al- 
lowed values of [3. To do this required asymptotic match- 
ing for the continuum solution (valid for n 1) to the 
"outer" solution of the difference equation with u> = 
(valid for n <C 1/w). This matching is imposed in the 
overlap region 1 <C n< 1/w. To study this limit, we let 
s — e* in Eq. ( 7.37| ), which is then dominated by s 3> 1. 
Thus 



f 



x(n,/3) ~ v/nw/2— / ds ( 
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J/3 



-i/3 



(7.38) 



This gives, 
X (n,M) 



2i(3 



-(l+M)/2 



(nui)' 



i/3/2 



(noj) 



1/3/2 



Comparison to the outer solution, Eq. (5.16) 



X (n,M)~n-( 1+M '>/ 2 c+n 
then gives, as in section V 



-^/ 2 +c^/ 2 



0k = Ttk/z. 

C3 = (-1) 



fc+1 C4 = C. 



(7.39) 



(7.40) 



(7.41) 
(7.42) 



We must next determine the constants c and a(/3). 
Normalization requires 



2\C k \ 2 \a(P)\ 2 Ip = 1, 



where 



In = 



2w(cosht + cosht') 



The next step is to rescale the parameter t — > t//3, t' — » r' / (3 to give 



2/3 4 



dtdt' sin t sin t'sech 2 



t-t' 



2u(cosht/f3 + cosht'/l3) 



dn , , „ , o 

— X(n,(3 2 - 
n 



/2u;(cosh t+cosh t') 



f2uj (cosh t //3+cosh t'/0) 



(7.43) 



(7.44) 



(7.45) 



(7.46) 



We are interested in small u>, with j3 = irk/\ ln(w)|. In 
this limit, 



exponential in Eq. ( 7.46 ) therefore becomes 



• - ( c 1 ' 3 + c pb-w-ss, (7.47) 
and an identical result with t — > t' . The factor in the 



Mcoshi + cosh<V-0{° 0<i,t'<7rfc/2 

L oo otherwise 

(7.48) 

Taking this limit therefore acts simply to restrict the lim- 
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its of integration, and we have 
1 



2/3 4 



rrk/2 ft-t' 

dtdt' sini sint'sech 2 [ — 

2/3 



(7.49) 



Since the sech is sharply peaked around zero in the (3 — > 
limit, we may effectively set t' « t in the second sine to 
obtain 



'/3 = 



1 

2/3 1 



tt/c/2 



g?£ sin i 



tt/c/2 



d£ 'sech 2 



(7.50) 



/o Jo \ 2/3 

The i' integral is clearly proportional to /3, and perform- 
ing these integrations exactly gives, in the [3 — ► limit, 

7rfc 

^ = 77^- (7-51) 



2/3 3 



k\Fh 



ci 



/2,Zy 



00 



Thus, by choosing 



a(/3) = /3 = 



hit 



we obtain the same constant c = z, 
hard-wall solutions in section V. 



-1/2 



(7.52) 



as found for the 



We are now in a position to calculate the Fermion 
Green's function. As in section V, to use the decom- 
position in Eq. ( [5.2| ), we must calculate matrix elements 
of single- Fermion operators betwee n the gro und and ex- 
cited states. In general, using Eqs. (4.13,5.6), 



(-l f(f)(n,M)x{n,-M) + {-l) k+1 (j>{n, -M) X (n, M) 
(-lfl(M) + (-l) k+1 I(-M) 



(7.53) 
(7.54) 



where the integral 



I(M) 



7rfc 



dn 



(f>(n,M)x(n,(3)n M/2 . (7.55) 



To evaluate this int egral , we let t — > e* in the integral 
representation, Eq. ( 7.29 ), giving 



(n, M) = ci 



UjM/2 



dte- V ^ cosht e Mt , (7.56) 



Inserting this and Eq. ( 7.36 ) above, the n-integration can 
be readily performed, giving 



I(M) = 



nkuj M / 2 

liOZ^j j ... 



,-Mti 



dt\dt' A 



cos /3£ 2 



coshii + coshi 2 



(7.57) 



The limits M, (3^0 can be safely taken in the numerator 
of the integral. The final result is 



I(M) = 



T 3 kLU ~(l~M)/2 
714 



(7.58) 



Putt ing this into Eq. ( |7.54| ) above and t hence into 
Eq. (5^2), one recovers the final result, Eqs. ( |5.3l| pT35 ) 
quoted in section V, with a different value, A = 7r e /4, for 
the nonuniversal constant. 



Deforming the contour to obtain a real integral gives 

cosh7ra f°° dw 
* Ji 



APPENDIX E 



In this appendix, we so lve the difference equation for 
J x eigenstates, Eq. (5.14), in the appropriate sector for 
the fermion Green's function. Consider the generating 
function, 



(7.59) 



n=0 



Multiplying Eq. (5.14) by w n and summing gives 



(1 — w 2 )——ip{w, a) — wif)(w, a) — —2i(nb(w, a). (7.60) 
dw 

This is easily solved by separation of variables. One finds 

jj(w, a) = (1 + w)- 1/2 - ia (l - w)~ 1/2+ia , (7.61) 

choosing ipo( a ) — 1 to fix the overall constant. This can 
be inverted using the contour integral 



dw i)(w, a) 
2ni w 1+n 



,1+n 



(w + l)- 1 / 2 -"^ - i)-i/2+» + (_!)«(«, + 1)-V2+'«( W _ iyi/2 



(7.62) 



(7.63) 
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For large n ^> 1, this integral is dominated by w w 1, 
and can b e easily evaluated to give the result quoted in 
Eq. (E35h. 
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